Kinetic glass behavior in a diffusive model 
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Three properties of the Edwards-Anderson model with mobile bonds are investigated which are 
characteristic of kinetic glasses. First is two-time relaxation in aged systems, where a significant 
difference is observed between spin and bond autocorrelation functions. The spin subsystem does 
not show two-time behavior, and the relaxation is stretched exponential. The bond subsystem shows 
two-time behavior, with the first relaxation nearly exponential and the second similar to the spin 
one. Second is the two-temperature behavior, which can be tuned by bond dilution through the full 
range reported in the literature. Third is the rigid-to-floppy transition, identified as a function of 
bond dilution. Simple Glauber Monte Carlo evolution without extraneous constraints reproduces 
the behavior of classical kinetic simulations, with the bond (spin) degree of freedom corresponding 
to configurational (orientational) disorder. 
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Kinetic glasses are long-lived nonequilibrium systems 
of considerable technological importance. The basic 
insight into their microscopic origin was provided by 
Kauzmanni there is a slowdown in configurational re- 
arrangement, caused by obstruction of kinetic motion. 
Despite a consensus in this view and advances in the- 
oretically enlightened phenomenology of known glass 
formers* 2 ^ there is at present no unifying view of vitrifi- 
cation. This is in remarkable contrast with its practical 
simplicity: put honey into the freezer or let egg whites 
dry, and a kinetic glass appears without effort. 

Out of equilibrium, subdominant microscopic correla- 
tions multiply, and finding the one responsible for macro- 
scopic structural arrest is even more difficult than doing 
so for equilibrium transitions. Simulations of abstract 
model systems have developed into an important tool to 
extract generic behavior in such a situation, although 
of course they do not substitute for direct physical in- 
sight into real cases^^^ While many models have been 
found which exhibit some kind of slowed-down response, 
their relationship with one another, and with physical 
reality, is still the subject of continuous investigation. 

Molecular dynamics (MD) simulations of two- 
component Lennard- Jones fluids are intuitively closest to 
real systems. 9 They have been compared with the predic- 
tions of mode-coupling theory , 10 ' 11 ! 12 a high-temperature 
microscopic approach, which depends on the explicit in- 
troduction of a three-particle scattering term. 

More abstract, and schematic, are constrained kinet- 
ics simulations.— These are lattice gas models in which 
updates are not dictated by pure coupling to the thermo- 
stat, but are supplemented with special rules, standing 
in for higher-order correlations, which otherwise cannot 
be included in a random walk. 

A separate class of models, sometimes used to un- 
derstand vitrification, are models of spin glasses^ It 
requires some abstraction, or leap of faith, to identify 



the former microscopic physical spin with a hypotheti- 
cal "mesoscopic" order parameter, representative of the 
slowed-down configuratonal rearrangement. Neverthe- 
less, particular quenched systems have remarkable paral- 
lels with other models of kinetic vitrification, the reasons 
for which have lately become better understoodi 15 ' 16 ' 17 ' 18 
Monte Carlo (MC) simulations play an important role 
in many of these investigations. In the time domain, they 
simplify joint probabilities as p(xi, ti; X2, £2! %3, £3) = 
p(xi,ti;x2,t2) for times t\ > t% > £3, since the underly- 
ing random walk has no memory. Here x% is the system 
configuration and the times are measured with a resolu- 
tion At shorter than the thermalization time, but long 
enough for all correlations, retaining information of the 
initial conditions, to die out, except the two-time ones. 
Such a random walker is a minimal model of fluctuations 
in a concrete thermostat. Since physical three-particle 
correlations are generally three-time correlations, in en- 
semble language one may say MC simulations take two- 
particle physical correlations as input and build higher 
correlations statistically as output, on scales coarser than 
At. 

In this work we address the issue as to how much needs 
to be said physically to obtain vitrification statistically. 
From this point of view, all of the approaches above suf- 
fer from some surfeit. Full Newtonian evolution does all 
the work physically, with no lower limit on At in princi- 
ple. Mode coupling has a three-particle kernel as input. 
Constrained kinetics has memory, since an update can 
overrule the thermostat. Conversely, to enforce the rules 
by the energy balance requires multibody forces, stronger 
than two-body. Random quenched disorder cannot be 
the input, since the comparatively well-organized vitre- 
ous disorde r 19 ' 20 is a consequence of evolution and by a 
Hamiltonian which has an ordered ground state. 

Our main result is that less needs to be said than sus- 
pected so far. We apply a natural annealment dynamics 
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to the Edwards- Anderson model for spin glasses,? 1 While 
the equilibrium state of the model was describe d 22 i 23 long 
ago, it turns out the approach to equilibrium is meaning- 
ful for kinetic glasses. When nearest-neighbor bond dif- 
fusion is added to spin flips, vitreous delay appears in the 
ordinary random walk. The difference with most other 
theoretical approaches is that there are two different de- 
grees of freedom in interaction, so the physical correlation 
responsible for vitrification appears to be a second-order 
off-diagonal term (mixed in conjugate fields). In contrast 
to other spin-glass models, it turns out that the bond, 
not spin, degree of freedom corresponds to the configu- 
rational one. Bond movement gives the model an unex- 
pected "off-lattice" character. 

We study the two-dimensional short-range Edwards- 
Anderson model,— as before , 24 i 25 at B = 0: 

H = -^J ij S i S j -M-B, (1) 

(i,3) 

where Si = ±1 and we shall use both a Gaussian distri- 
bution of Jij with half-width J and a bimodal ± J distri- 
bution. We investigate a two-dimensional lattice of size 
500 x 500 with Glauber dynamics. Between each two spin 
trials (flips) is a bond trial: a positive and negative bond 
impinging on the same site, and chosen at random, are 
allowed to exchange places, subject to the same criterion 
as the spin trials. The concentration paf of antiferro- 
magnetic (AF) bonds is taken to be 50%. 

Bond updates quickly anneal the sample, and after 
a transient of ^100 updates per site it enters a long- 
lived state with glasslike dynamics,— and a second, 
much longer relaxation time^ The "glassiness" of the 
metastable spin state is not topological: the underlying 
bond distribution has low frustration at any instant in 
time, so much so that it can be mapped onto a disordered 
ferromagnetf 2 ^ with a finite transition temperature^ 
Bond diffusion neverthless prevents the spins from set- 
tling into any given ordered state, leading to a decay of 
correlations even below the hidden phase transition. 

When bonds move, the model theoretically evolves to- 
wards an annealed equilibrium! 22 ' 23 However, as long as 
bond updates are kept local, their diffusion is not very 
efficient in finding the optimal configurations, despite ef- 
ficient annealment in energy. We do not establish directly 
that the configurations are suboptimal, but when bond 
hops of arbitrary range are allowed, the configurations 
obtained are quite different than those found here, hav- 
ing a strong tendency to dropletlike phase separation. 2 - 
Although equilibration is delayed from the point of view 
of correlations, evolution passes only through a small sub- 
class of energetically favorable states, though bond up- 
dates imply that not only the equilibrium spin manifold 
is being sampled. This is similar to the situation in real 
kinetic glasses, while in spin glasses the average config- 
uration energy is much higher^ We find that isotropic 
spatial correlation functions fall off rapidly within a few 
lattice spacings, in accordance with recent thinking. 29 




FIG. 1: Spin and bond autocorrelation functions in an aged 
sample. In each pair of curves, the bond autocorrelation is the 
one with the plateau at intermediate times. The full, dashed, 
dotted, and dot-dashed curves correspond to temperatures 
T/J = 0.05, 0.15, 0.30, and 0.50, respectively. Each Monte 
Carlo step (MCS) consists of a spin and a bond update per 
site. 

In Fig. [TJ we show the autocorrelation functions of two 
copies of an aged sample, for spins given by 

i v 

C{t) = -Y,{S i {t )S i {t Q +t)), (2) 

i— 1 

and similarly for bonds, for the Gaussian distribution. 
"Aged" means to long enough for no transients to sur- 
vive after quenching, so the system is in the long-lived 
metastable state. The spin and bond autocorrelation 
functions are different. The spin autocorrelation decays 
as a stretched exponential exp [— (i/r)' 3 ], which is the 
usual behavior in disordered spin models. The bond au- 
tocorrelation shows obvious two-time behavior, with a 
plateau separating the fast, closer to exponential decay 
(f3 w 0.9), from a slow, stretched exponential component 
(f3 f» 0.8). The plateau is characteristic of configura- 
tional relaxation in kinetic glasses and observed in clas- 
sical simulations of binary systems. 9 The fast component 
"shakes out" those configurations which may be relaxed 
locally, while the slow component refers to "locked-in" 
ones, whose relaxation is impeded by intermediate-range 
correlations. We recall^ that the orientational degree of 
freedom is generally not affected by the glass transition, 
being slowed down already above T g . As the tempera- 
ture is lowered, the spin and bond responses decouple 
ever more, which is visible in the figure as an extended 
plateau in the bond response. This is also analogous to 
real systems, where there are more translations relative 
to rotations the deeper one goes into the glassy regime. 3 
If spin updates are discontinued, the bond response loses 
the second relaxation component, and the plateau then 
extends to the longest times investigated. The reason 
is that fixed annealed spins create barriers of satisfied 
bonds, which mobile bonds find difficult to cross, so the 
system is broken into uncommunicating regions. When 
spin updates are allowed, these regions themselves evolve 




FIG. 2: Violations of the fluctuation-dissipation theorem for 
an (a) Gaussian and (b) bimodal distribution, (a) Bottom to 
top: T/J = 0.05, 0.01, 0.015, 0.2, 0.25, 0.3, 0.4, 0.5, 0.6, respec- 
tively, (b) Top to bottom: p = 100%, 90%, 80%, 70%, 60%, 
40%, 30%, 20%, 10%, respectively, at T/J = 0.2. 



slowly, accounting for the second relaxation time. Rare 
spin updates at low temperature correspond to activated 
hops in constrained kinetics, absent from mode coupling, 
which similarly retains the plateau ad infinitum. 

Next, we investigate how the interplay of the two de- 
grees of freedom affects departures from the fluctuation- 
dissipation theorem (FDT), expected in the same 
context.-^ Here we plot C(t) for spins [Eq. (J2J)] against 
the spin susceptibility normalized by temperature: 



, s T . AM 
T x (t) = — lim - — 
Ay ' V ab^q AB 



(3) 



t +t 



for AB(to + t) a step function at t = 0. Although such 
violations are well established for models with quenched 
disorder, it is not clear that annealing caused by bond 
diffusion will not quickly lead to equilibration. In fact, 
it does not, and for the same Gaussian bond distribution 
as above, we obtain in Fig. [2] the "many-temperature" 
violation curve, known from previous investigations of 
spin models and observed in some experiments. 30 Since 
the model is only coupled to a simple thermostat, the 
violation is evidently generated internally. 

Even more interesting is the case of the ± J bond dis- 
tribution. We introduce bond dilution through a fixed 
proportion p$ of bonds with J = 0. For po = 1, an or- 
dinary paramagnet, the FDT is evidently observed. As 
the dilution is reduced there are again violations, but 
instead of many temperatures, the "two-temperature" 
shape appears, familiar from MD simulations of two- 
sphere models. 31 ' However, in the ±J case there is no 



FIG. 3: Rigid-to-floppy transition. The fraction of indifferent 
spins is shown as a function of bond dilution. The minimum 
corresponds to the least-strained system. The inset shows the 
minimum as a function of temperature. 



plateau in the bond response, indicating that a fine en- 
ergy scale is needed for the creation of inherent struc- 
tures. 

The rigid-to-fioppy transition has emerged as an im- 
portant paradigm in a wide range of glass formers4^ 
The basic ide a 19 ' 20 is that a glassy system can be ei- 
ther overdetermined (rigid) or underdetermincd (floppy), 
with some ideal "unstrained-but-stiff" configuration in 
the middle. Figure [3] shows this model" s version of the 
rigid-to-floppy transition. It gives the fraction p of spins 
in a zero local field ( "indifferent" ) as a function of bond 
dilution pq. As po — > 1 (paramagnet), all spins become 
indifferent: p(l) = 1 is the extreme floppy limit. How- 
ever, the curves p(po) also show a sharp minimum at a fi- 
nite pq. This is not difficult to understand: an indifferent 
spin in an undiluted system occurs wherever two positive 
and two negative bonds impinge on the same site. It can 
be oriented by putting a single zero bond in place of one 
of them. Hence the number of indifferent spins initially 
decreases with dilution. The minimum corresponds to 
the least-strained state in the sense above: to the left, p 
increases because there are too many constraints and to 
the right, because there are too few. The trend of the 
minimum with temperature is also reasonable; one ex- 
pects that the lower the temperature of the quench, the 
more zero bonds are needed to effectively unstrain the 
system. 

Experimental evidence suggests that vitrification is 
a two-stage process. The first stage is dominated by 
quantum mechanics — i.e. chemical bonding — which de- 
termines what is "locally favorable." Once the first stage 
has given rise to entities at the scale of ~ 10-100 nm, 
their further evolution may be imagined in classical or ge- 
ometrical terms. Most investigations study this second 
stage of vitrification, including the present work. The 
paradigm of a macroscopic internally generated poten- 
tial, borrowed from phase transitions, is evidently sensi- 
ble, since something is certainly precluding window glass 
from free flow. However, the universally observed liquid 
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structure factor— means this potential varies in a way 
irreconcilable with a divergent length scale. Thus the 
notion of an internally generated field is more general 
than that of a phase transition. In addition, macroscopic 
approaches f^i^i based on dynamic scaling ideas, indicate 
that the time scales of the internal field are finite as well. 

Given that classical simulations do not at present cor- 
respond to any well-defined physical objects — indeed the 
"configurational degree of freedom" is as intuitive a con- 
cept today as when its importance was first noted^ — 
it is reasonable to ask whether their success in repro- 
ducing particular observed aspects of glassy dynamics 
comes from the side of the thermostat, rather than be- 
ing a higher-order effect of precise mechanical evolution. 
From this point of view, ordinary MC simulations are a 
fundamental tool. They build up high-order correlations 
without microscopic trajectories, enabling one to study 
how close to equilibrium glasses really are. 

The present work obtains a vitreous response with a 
minimal disturbance of the annealed spin equilibrium and 
the minimal physical input to date. All that seems re- 
quired is dynamical evolution of two different degrees of 
freedom with a local interaction. Hence vitrification in 
the model is due to a second-order off-diagonal correla- 
tion, rather than a higher-order diagonal one, typical of 
models with one kind of degree of freedom. It raises the 
possibility that the same could be true in real systems, 
such as large molecules, the interdependence of trans- 
lation and rotation (both kinetic) in a log-jam coming 
readily to mind. The idea of configurational rearrange- 



ment agrees with the interpretation of off-diagonal cor- 
relations as conditional probabilities: a particle cannot 
move unless another one does. However, the definite 
identification of the correlation involved depends on find- 
ing a natural conjugate field, a nontrivial task, even in 
a simple model such as ours. Reproducing MD by pure 
MC simulations^ is very similar to the present result; 
we have removed from it one more physical prop, the re- 
alistic potential, so it only seems important that there 
be two kinds of little balls. This fits well with binary 
(and ternary) mixtures being easier to vitrify than pure 
substances 

Another issue is that the model has a hidden ferro- 
magnetic transition . 26 ! 27 It is analogous to the avoided 
crystallization in vitreous liquids. The bond autocorre- 
lation implies that discontinuing spin updates triggers a 
true transition, as is known to happen in mode coupling. 

To conclude, we have given a concrete example of vit- 
rification in a simple model, where it is by default due 
to a second-order off-diagonal correlation. The appear- 
ance of vitreous slowdown in the bond response of the 
Edwards- Anderson model when bond diffusion is intro- 
duced points to its kinetic origin. We believe this type of 
correlation may be responsible for vitrification in at least 
some real glasses, where it would provide the minimal 
physical content of the configurational degree of freedom. 
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ernment under Project No. 0119256. 
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